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ABSTRACT 

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty 
torus around the central source, giving rise to Seyfert 1 line spectrum for pole-on 
viewing and Seyfert 2 characteristics in edge-on sources. Although the observed 
IR is in broad agreement with this scheme, the behavior of the 10 /im silicate 
feature and the width of the far-IR emission peak remained serious problems 
in all previous modeling efforts. We show that these problems find a natural 
explanation if the dust is contained in ~ 5-10 clouds along radial rays through 
the torus. The spectral energy distributions (SED) of both type 1 and type 2 
sources are properly reproduced from different viewpoints of the same object if 
the optical depth at visual of each cloud obeys r v > 60 and the clouds' mean 
free path increases roughly in proportion to radial distance. 

Subject headings: dust: extinction — galaxies: active — galaxies: nuclei — galaxies: 
Seyfert — quasars: general — radiative transfer 

1. INTRODUCTION 



Although there is a bewildering array of AGN classes, a unified scheme has been emerg- 
ing steadily (e.g. Antonucci 1993, 2002; Wills 1999). The nuclear activity is powered by a 
supermassive (~10 6 -10 9 M & ) black hole and its accretion disk, which extends to ~ 1 pc. 
This central engine is surrounded by a dusty toroidal structure, extending to > 100 pc. 
Much of the observed diversity is simply the result of viewing this axisymmetric geometry 
from different angles. The torus provides anisotropic obscuration of the central region so 
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that sources viewed face-on are recognized as Seyfert 1 galaxies, those observed edge-on are 
Seyfert 2. The primary evidence for the torus comes from spectropolarimetric observations 
of type 2 sources, which reveal hidden type 1 emission via reflection off material situated 
above the torus opening. While compelling, this evidence is only indirect in that it involves 
obscuration, not direct emission by the torus itself. 

An obscuring dusty torus should reradiate in IR the fraction of nuclear luminosity it 
absorbs, providing direct evidence for its existence. Indeed, the continua from most AGN 
show significant IR emission. Silicates, whose presence is revealed by a spectral feature at 10 
/iin, are a major constituent of astronomical dust. Seyfert 2 sources display this feature in 
absorption — as expected for edge-on viewing of an optically thick torus. However, contrary 
to expectations for face-on viewing, quasars and Seyfert 1 galaxies do not show this feature in 
emission, (Roche et al 1991; recent claims by Clavel et al 2000 notwithstanding). This poses 
a serious problem to models of the dusty torus emission, which must suppress the silicate 
feature in type 1 objects while producing it in type 2. Laor and Draine (1993) studied the 
effects of altering the dust properties and conclude that suppression of the 10 /im feature 
requires significant silicate depletion or very large (up to 10 /im) grains (see also Maiolino et 
al 2001). This could explain type 1 SED but would require a different explanation for type 
2 objects with their prominent 10 /im absorption feature. 

Pier & Krolik (1992, 1993) were the first to explore the effects of toroidal geometry on 
dust radiative transfer. They note that the AGN dust must be in clumps to protect the grains, 
but because of the difficulties in modeling a clumpy distribution they approximate it with 
a uniform one instead. They also neglect scattering. In spite of these approximations, their 
results are encouraging. The directional dependence of their model radiation reproduces the 
gross features of observed SED, indicating that the toroidal geometry captures the essence of 
IR observations. Two major problems remain: the observed emission has a far-IR peak much 
broader than the models can produce, and the silicate emission feature is suppressed in face- 
on viewing only in a narrow, finely-tuned range of the model parameters. Rowan-Robinson 
(1995) notes that both problems could be alleviated by dumpiness. He argues that a spherical 
shell is an adequate proxy for a clump, and produces type 1 IR emission from a superposition 
of such concentric shells. But the shell/clump analogy is fundamentally problematic: Since 
the illuminated and dark faces of a clump produce widely different spectra, the observed SED 
varies with the angle between radiation source, clump and observer (see fig. 1). In contrast, 
outside an optically thick spherical shell the SED is always the same irrespective of observer 
location because only the dark side is ever visible (inside the shell, only the illuminated face 
is visible). Concentric spherical shells cannot substitute for clumps. 

We performed a more adequate calculation of IR emission from an obscuring clumpy 
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torus, combining realistic modeling of the emission from externally illuminated clumps with 
proper handling of cloud shadowing. The latter is an essential ingredient since the number 
of clouds along radial rays through the torus must be sufficiently large to ensure that x-rays 
from the central engine are virtually always attenuated in type 2 objects (Guainazzi et al 
2001). Here we report the results of our effort, which show that the SED of both type 1 
and type 2 objects are properly reproduced at different viewing angles of the same clumpy 
torus — in agreement with unified schemes. In a forthcoming publication we will provide a 
more detailed account of our calculations and the results. 



Consider a medium where the dust is in clouds. For simplicity, each cloud has the same 
optical depth t\. Along a given path, the mean number of clouds encountered in segment 
ds is dj\f(s) = ds/£(s), where £ = (n c A c ) _1 is the mean free path between clouds (A c is 
the cloud area perpendicular to the path and n c is the number density of clouds). As long 
as £ ^> R c , where R c is the cloud radius, each cloud can be considered a point source of 
intensity S c \, and the intensity generated along the segment is S c \dJ\f. Denote by P esc the 
probability that this radiation propagate along the rest of the path without absorption by 
any other cloud. Natta & Panagia (1984) show that Poisson statistics for the distribution of 
clumps yields 



of segment ds to the emerging intensity is simply P esc S c \dJ\f and the flux from the cloud 
distribution at distance D is 



where dA is the surface area element perpendicular to the line of sight. Given the geometry, 
the flux can be calculated from this expression once S c \ is known. 

Clouds are heated by radiation from both the AGN and all other clouds. Consider first 
the direct heating by the AGN (figure 1). Since our interest is in optically thick clouds only, 
the dust temperature is much higher on the illuminated face than any other part of the 
surface. Therefore the cloud emission varies strongly with direction, and the corresponding 
source function S^ x (r, a) depends on both distance r and position angle a. The clump shape 
is of course arbitrary, and we have constructed "synthetic clumps" by averaging the emission 
from an illuminated slab over all possible slab orientations. This procedure utilizes an exact 
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solution of radiative transfer for external illumination and also reproduces the a-dependence 
of the observed fraction of illuminated area on the surface of a spherical-like object. 

We conducted detailed calculations with the code DUSTY (Ivezic, Nenkova & Elitzur 
1999). The code preforms an exact solution of the slab radiative transfer problem including 
dust absorption, emission and scattering. The optical depth across the slab thickness at 
wavelength A is t x = q\Ty where ry is the optical depth at visual and q\ is the proper 
efficiency factor of standard interstellar dust. With AGN luminosity L 12 = L/1O 12 L , the 
illuminating bolometric flux is F e = L/Anr 2 . Its spectral shape, typical for AGN, is \f\ 
= constant in the wavelength range 0.01-0.1 jum and oc A -0 5 for 0.1-100 jum (cf Pier & 
Krolik 92, Laor & Draine 93). For any slab orientation, the solution of the radiative transfer 
problem determines the temperature run in the slab and the intensity it emits toward any 
direction. Therefore the dust temperature T n on the illuminated face of a slab normal to the 
radius vector can replace the external flux F e as a specifier of location. 

The top panel of figure 2 shows typical results for Sf x . The clump spectrum was con- 
structed by averaging over all slab orientations the solutions for slabs with ry = 100 at radial 
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distances where T n = 800 K (corresponding to r = 4/^ pc). The 10 fim silicate feature 
is seen in absorption at small a, switching to increasingly more pronounced emission as a 
increases and a larger fraction of the clump illuminated face becomes visible. 

An exact calculation of the effects of diffuse radiation from all clouds is a formidable 
task, akin to a full solution of the radiative transfer problem in which individual dust particles 
themselves are dusty clouds. However, compared with the AGN radiation, heating by the 
diffuse radiation is highly inefficient because it involves long wavelengths, and thus can be 
neglected. Emission from clouds with direct view of the AGN is well described by S^ x . 
But clouds whose line-of-sight to the AGN is blocked by another cloud will be heated only 
indirectly by the diffuse radiation. We approximate the diffuse radiation field at r by angle- 
averaging over a the emission Sf x (r,a) from clouds directly heated by the AGN. In this 
radiation field we imbed a sphere with constant density profile and solve for its temperature 
and emission. Our approximation for the source function of indirectly illuminated clouds is 
Sl x = F/Q, where F is the flux and Q the solid angle of such a sphere at a large distance. 
The bottom panel of figure 2 shows a sample of S 1 ^, together with the corresponding source 
functions for direct heating. 

At distance r the probability for unhindered view of the AGN is p(r) = e~ A/ ' (r - ) , where 
M(r) = f dr/£ is the mean number of clouds to the AGN. Our approximation for the clump 
source function is 

S c x=pS* x + (l-p)Si x . (3) 
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The steps we outlined can be iterated to become an exact solution scheme for radiative trans- 
fer in clumpy medium. The small magnitude at short wavelengths, which control heating, 
of Sl x / S^x evident in figure 2 indicates that this process should converge rapidly. Eq. 3 is 
equivalent to its first two steps. 



3. MODEL RESULTS 

We model the AGN obscuring region as a toroidal distribution of dusty clouds, shown 
in figure 1 . Instead of the inner radius R\ we specify as input T n ; = T n (R{), the surface 
temperature of a normally illuminated slab at R\. In all calculations we set T ni = 1400 K, 
the temperature below which both silicate and graphite grains exist; this choice implies R\ 

1/2 

= 1.2 L 12 pc for a Ty = 100 cloud. The geometry requires two additional input parameters: 
The torus thickness Y = R Q /Ri, equivalent to the lowest T n in the torus, and 9, the conical 
opening half-angle. Assuming power law variation for the mean free path £ cx r q , the cloud 
distribution is described by the two input parameters q and J\f T = drji, the number of 
clouds along a radial ray through the torus. The final free parameter is Ty, the optical depth 
of each cloud. The grain efficiency factor qx and the AGN spectral shape fx are set to their 
standard values. The observed flux is determined from 

where F x is calculated from eq. 2 and txT from eq. 1 with Af(s) = A/t for type 2 sources 
and zero for type 1. Extracting the overall bolometric flux Fb = L/AnD 2 , the SED is 
controlled by T ni and five dimensionless free parameters: Y and specify the geometry, Ty 
the individual clouds and A/t and q the cloud distribution. 

We performed extensive calculations and find a large range of the free parameters that 
yield satisfactory agreement with observations. A detailed report of our results will be 
provided elsewhere, here we show in figure 3 a representative model. Our aim is to reproduce 
the typical SED of the type 1 and type 2 families, not to fit in detail any particular object. 
Data points for type 1 sources show average spectra for radio-quiet quasars (Sanders et al 
1989, Elvis et al 1994) and Seyfert 1 galaxies (Granato & Danese 1994). Because of the 
high obscuration of the AGN in type 2 sources, fluxes for their nuclear regions properly 
extracted out of the contributions of the host galaxy and starburst regions are scarce. We 
show data for reliable measurements of four individual objects. The data were scaled for a 
rough match of the model results without attempting a best fit. Detailed fitting of type 2 
objects, when warranted, will determine their bolometric flux Fb, a quantity inaccessible for 
direct measurement because the bulk of the flux is emitted away from Earth. 
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Our models reproduce the broad IR bump extending to ~ 100 /im, as observed. Fur- 
thermore, the 10 /iiii feature appears in absorption in equatorial viewing but is smeared out 
in axial viewing in spite of its prominence in emission from directly-heated individual clouds 
(fig. 2). The reason is that most of these clouds are blocked from view even along the axis. 
We find satisfactory results for conical opening G = 45 ± 15°. The torus thickness can vary 
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in the range Y = 50-250, i.e., R Q /L 12 = 60-300 pc. At smaller Y the IR bump is not 
broad enough, larger Y produce too much 20-100/zm emission. The number of clouds is A/t 
= 5-10. When Af T < 5 the 10/xm feature appears in emission in pole-on spectra, J\f T > 10 
shifts the spectrum to longer wavelengths resulting in too little 1-10/mi emission and too 
much 20-100/xm. The power index should be q = 1 ± 0.5. At larger q the far-IR bump is 
not high enough, as is evident in figure 3, at smaller q it is too high. 

Significantly, the only constraint on the optical depth of individual clouds is r v > 60. 
The results vary only slightly when r v increases from 60 to 100, and hardly at all during 
further increase. The reason is simple. The dependence on Ty arises from the probability for 
photon escape and the cloud source function. From eq. 1, P csc = whenever t\ 3> 1, and 
at ry > 60 this condition is met at all relevant wavelengths. Similarly, because the clouds 
are heated externally, only their surface is heated significantly when Ty is large. Increasing 
Ty further only adds cool material, thus S c \ saturates for all relevant A, similar to standard 
black-body emission. Extending our calculations all the way to r v = 500, we have verified 
that increasing Ty indeed has no effect on the model results. 

4. DISCUSSION 

The cloud distribution is described by q and A/"t, individual clouds by the optical depth 
ry. No other cloud property was specified. The cloud size enters only indirectly through the 
underlying assumption R c <C I. If V c ~ A C R C is the volume of a single cloud, the volume 
filling factor of the cloud population is = n c V c = R c /£. Our calculations apply to small 
filling factors, <C 1, the Pier & Krolik model involves the opposite limit 0=1. The mean 
free path scale can be determined from £(Ri)/Ri = (1/A/r) y~ q dy. The model with q = 

1 /2 

1, A/t = 5 and Y = 100 has £{R\) = 0.9/?; = 1.1 L 12 pc. A reasonable realization of this 
model, though not unique, is a constant = 0.1 throughout the torus so that R c = 0.09r. 
When ry = 100 and L 12 = 1, the clouds vary from R c ~ 0.1 pc with gas density ~ 3 10 5 
cm -3 in the torus inner regions to R c ~ 10 pc with density ~ 3 10 3 cm -3 at the outer edge. 
The torus could contain additional clouds with smaller ry without significantly affecting the 
SED. 

The properties we deduce for the torus agree with independent estimates. Based on 88 
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Seyfert galaxies Schmitt et al (2001) conclude that = 48°, at the center of the range we find. 
Our minimal cloud number At = 5 implies a probability < e~ 5 = 7 10~ 3 for unattenuated 
view of the AGN in type 2 sources, in agreement with the findings of Guainazzi et al (2001). 
From x-ray measurements of 73 Seyfert 2 galaxies Bassani et al (1999) find a large variation 
in column density, with a mean of 3 1 23 cm~ 2 and > 10 24 cm~ 2 for as many as a third of 
the sources. The mean is comparable with the column density of our torus for 5 clouds of 
r v > 60 each and standard dust-to-gas ratio. Furthermore, the similarity of SED among 
type 2 sources in spite of the large x-rays column variation is a natural consequence of our 
model: Ty is bounded only from below, the SED remains the same at all A > 1/iin for 
arbitrary increase in overall column. But x-ray absorption does vary with the latter because 
the optical depth for Thompson scattering is only ~ 10~ 2 rv. Sources with small columns 
may show up in x-ray absorption while selectively excluded from IR observations because of 
their weak emission. 

Our results add strong support for unification schemes of AGN. In accordance with 
such schemes, the IR emission from both type 1 and type 2 sources is reproduced at different 
viewing of the same geometry when proper account is taken of the torus dumpiness. 

We have greatly benefitted from discussions with many colleagues, especially J. Conway, 
A. Laor, N. Levenson and H. Netzer. Support by NASA and NSF is gratefully acknowledged. 
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Fig. 1. — Model geometry. Top: Positions of the AGN, clump and observer. As the position 
angle a varies at fixed distance r, the visible fraction of the clump's illuminated area changes 
and with it the observed radiation. Bottom: The clumpy torus. 
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Fig. 2. — Emission from ry = 100 clumps normalized to the AGN local bolo metric flux 
F e = L/Aur 2 . Top: The source S^ x for directly heated clouds at distance where T n = 800 K 

1/2 

(r = AL 12 pc for this ry). The position angle a is shown at 22.5° intervals. Bottom: Emission 

of directly (S^ x , full lines) and indirectly (S^ x , dashed lines) heated clouds at distances where 
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T n = 1400, 500 and 200K, or r/L^ = 1, 10, 100 pc, respectively. The source function for 
direct heating is shown at a = 90. 
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Fig. 3. — Modeling and observations of type 1 (top panel) and type 2 (bottom) sources. 
Lines are model results for the axial and equatorial views of the clumpy torus shown in 
figure 1, with F B = L/AttD 2 . Radial rays through the torus have 5 clouds on average, each 
with Ty = 100. The clouds' mean free path varies as r q , with the indicated q. The torus 
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inner radius corresponds to temperature T ni = 1400 K (i?i = 1.2 L 12 pc) and it has R D = 
100i?i and G = 45°. Data for type 1 are composite spectra from the indicated compilations, 
type 2 are for individual objects. Each data set was scaled for rough matching with the 
models. 



